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The critical properties ol an infinitely long Ising strip with finite width L joined periodically 
' or antiperiodically are investigated by analyzing the distribution of partition function zeros. For 

("S , periodic boundary condition, the the leading finite-size scaling of partition function zeros and its 

corrections are given. For antiperiodic boundary condition, the critical point of 2D Ising transition 
is one of the loci of the zeros, and the associated non-analyticity is identified as a first-order phase 
, transition. The exact amount of the latent heat released by the transition is 4/L. 

PACS. 05.50.+q - Lattice theory and statistics. 
' pV ' PACS. 75.10.Hk - Classical spin models. 

, PACS. 75.70.Cn - Magnetic properties of interfaces. 
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The geometry of the system has been known to play a crucial role in many aspects of critical behavior: The 
+L , way of taking the thermodynamic limit may affect the occurrence of a phase transition, and critical exponents 
depend on spatial dimensionality non-trivially. These result in more stringent test on the universality by including 
universal critical amplitudes and amplitude relations obtained from the finite-size effectsP,0- Boundary condition is 
a important factor in determining the finite-size scalings. In particular, aperiodic boundary conditions may introduce 
interfaces to the systems and lead to a profound change in the scaling behavior 0, 0]- In this letter, we report the 
appearance of a new first-order transition induced by the interface of a 2D semi-infinite Ising system. We consider the 
' simplest system, an infinitely long Ising cylinder defined on a square L x x L y lattice with infinite L x — > oo and finite 
L y = L. To have the system with cylindrical geometry, we impose periodic or antiperiodic condition on the infinitely 
extended boundary rows, (m, 1) and (m,L) with — oo < m < oo. The condition on the couplings of the boundary 
spin variables cr mj i along the y-axis is o~ m ,L+i = &m,i for periodically joined circumference (PJC, hereafter), and 
the condition becomes er m ,L+i = -fm,i for antiperiodically joined circumference (AJC, hereafter). The free energy 
' of the system with AJC contains the additional contribution from the interfacial tension in comparison with that of 
. the system with PJC0,|(|. To understand the differences in thermodynamic properties between the two boundary 
y—i ' conditions, we analyze the distribution and the density of partition function zeros. Our results reveal some intriguing 
C*~) , properties for the system with AJC of finite L. 

We already have a complete picture about the partition function zeros of the zero-field square Ising model in the 
thermodynamic Iimit0,l3- Fisher first showed that the partition function zeros, referred as Fisher zeros, all lie on the 
unit circle, \z\ = 1, in the complex plane of z = sinh (2/3) with the zero on the positive real axis, z c = 1, corresponding 
to the critical point of Ising transition^. Then, Lu and Wu obtained the density associated with the zeros on the 
•p-; , unit circle analytically^. Based on these results, we study the finite-size scalings of Fisher zeros for infinitely long 
Ising cylinders. For the system with PJC, portions of the zeros, including the critical point, are always absent from 
the unit circle for finite L, and the scaling behavior is given in the exact form. However, for the system with AJC, 
regardless the L value, the zeros distribution always contain the critical point. Then, by analyzing the cumulative 
distribution of Fisher zeros around the critical point we show that the singular property of the critical point can be 
classified into two types: One is the first-order transition for a finite circumference, and the other is the second-order 
transition for an infinite circumference. The latter corresponds to the Ising ferromagnetic phase transition, and the 
former is a new phase caused by the fluctuation of interface from AJC. Thus, in taking the thermodynamic limit 
L — > oo from an infinitely long Ising cylinder of finite L, we encounter different situations for two different boundary 
conditions: (i) For the case of PJC, only the second-order phase transition appears in the limit L — > oo, and there is 
no phase transition in the intermediate stage < L < oo. (ii) For the case of AJC, we first have the first-order phase 
transition which becomes weaker as L increases and eventually reduces to a second-order phase transition in the limit 
L — > oo. This picture is further verified by measuring the latent heat per site released in the first-order transition as 
4/L. 

The partition function of a plane square L x x L y Ising lattice with periodic or aperiodic boundary conditions can 
be expressed as the sum of four terms 0,0, but the sum reduces to one term in the limit of L x (or L y )—y oo. In this 
limit, the system becomes an infinitely long Ising strip with finite width L, and the corresponding free energy of the 
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for the case of isotropic ferromagnetic couplings|5j,|6|]. Here the possible p- values in the sum are half integers ranging 
from 1/2 to L — 1/2 for PJC and integers ranging from to L — 1 for AJC, and the function $p(</>) is given as 

27TD 

$ p (0) = cos0 + cos — — . (2) 
Then, the loci of the Fisher zeros can be obtained as the union of the solutions of the condition 

z+ --%{(j>) = (3) 

z 

with < 4> < 2n for all allowed p-values. The solutions of Eq. (3) are expressed as 

z = exp (±i0) , (4) 
for < 9 < tv, where 9 are the elements belonging to the union of the sets of 9 p for all allowed p-values with 

cos p ^^M. (5) 

For the function $ p ((f>) of Eq. (2), the value of 9 P for a given p-mode is in the range between P) mi n = 
cos -1 [cos 2 (np/L)] and p , ma x = 7r — cos -1 [sin 2 (np/L)]. An example of the curve of 9 P (</>) obtained from Eq. 
(5) for various p-modes with L — 6 is shown in Fig. 1. 

For PJC, by including all allowed p-modes we obtain the range of 9 as 9 m - m = 9 p=1 / 2y m in an d 9 max = tt — 
cos -1 [cos 2 (7r/2X)] for even L, and 9 min — p= i/2, m in an d ^max = ?r for odd L. Thus, as shown in Fig. 2(a) for 
L = 5 and 2(c) for L = 6, the zeros distribution can be summarized as the followings: (i) The zeros are absent from 
two arcs of the unit circle for even L and from an arc for odd L, and the zeros do not contain the critical point for 
both even and odd finite L. (ii) By denoting the nearest zero to the critical point as z\ (L), we obtain the opening 
angle between z% (L) and the critical point as (9 m i n . To obtain the finite-size scaling behavior of z% (L), we expand 
sin 9 m i n and cos m i n as a power series of 1/L to obtain 

^ r /IT \2fe+l 

Im^(L) = ^C| fe+1 (— ) , (6) 



fc=0 



with the coefficients C{ = y/2, C| = -5^2/12, C{ = 49%/2/480, and etc.; and 



\Rez 1 (L)-z c \=Y / C«(^-y\ (7) 



fe=i 



with the coefficients C 2 = 1, = —1/3, Cg = 2/45, and etc.. From the finite-size scaling theory of the zeros 
ditributions, the leading finite-size scaling behavior of the imaginary part of a Fisher zero labelled by j is given as 

Imz^Lj-r 1 /", (8) 

where v is the correlation-length exponential; an d the leading scaling behavior of the real part of the lowest zero 
(j = 1) can be written as 

|Rezi(£)-* c | ~L" A — , (9) 

where A zero is another critical exponent which is closely related to the shift exponent A characterizing the shift of the 
specific-heat peak from the critical point^J- Then, from Eqs. (6) and (7) we have v — 1 and \ zer o = 2. The two 
values are the same as the results obtained by Janke and Kenna for a finite rectangular lattice with the Brascamp-Kunz 
boundary conditional but different from those obtained by other boundary conditions [l2T|. 

For AJC, we obtain the range of 9 as < 9 < 6> max , with 9 max — it for even L and 9 max = ir — cos -1 [cos 2 (7r/2£)l 
for odd L. Thus, as shown in Fig. 2(b) for L = 5 and 2(d) for L = 6, the zeros distribution is absent from an arc of 
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the unit circle for odd L and fills the unit circle completely for even L. Moreover, the critical point is one of the loci 
of the zeros for both even and odd L. Then, a question naturally arises: What is the non-analytic property associated 
with the zero at the critical point? This can not be the second-order phase transition which occurs only in the limit 
L — > oo. To determine the nature of the non-analyticity, we study the functional form of the cumulative distribution 
of the zeros [l3] . 

The cumulative distribution of zeros, Gl (9), is defined as the total number of Fisher zeros in the interval [0, 8} of 
the unit circle. By expressing in terms of the zeros density ql (0), we have 



Gl{9)= f 9 l 
Jo 



(x)dx. (10) 



Here we always take the normalization as Gl (t) = 1/2 for any width L. It has been shown [Til EUthat the zeros 
density <?oo (9) near the critical point, 9 C — 0, behaves as (9) a29 la in a second-order phase transition for 
which the specific heat diverges as c ~ \T — T c \ a . Thence, the corresponding cumulative distribution of zeros is 
Goo (9) ~ 62# 2 ~°- In 2D Ising transition for which we have a — 0, the density goo (9) was shown to vanish linearly as 
5oo (9) « '\0\ /2ir near the critical point 8], and the corresponding Goo (9) obtained by the integration of g^ (9) given 
by Ref. Q are shown as a solid line in Fig. 3. For the first-order transition, the asymptotic behavior of <7oo (9) to the 
critical point was shown to take the form, goo{9) ~ 5oo(0) + a\6 n + with goo(0) proportional to the latent heat|l6|. 
Correspondingly, the behavior of Goo (9) is 

G oo (9)^g oo (0)9 + b 1 9 n+1 + .... (11) 

For the system of an infinitely long Ising cylinder, the value of Gl (9) is proportional to the sum of the lengths of 
the curves 9 P (cj>), as shown in Fig. 1, over all possible p- values up to a given 9 value with < cj) < 2ir. Because of 
the symmetry, 9 P (</>) = 9 P (—0), we may restrict the (^-values to the interval [0, 7r]. Then, by using the normalization 
condition Gl (tt) = 1/2 we write 

G *(*) = -r-E ( 12 ) 

<lltot — 

p 

where l p (9) is the length of the curve 9 P (<fi) up to a given 9 value, 




and Itot is the sum of the lengths of all the curves, 

hot = J]] lp (#p,max) • (14) 

p 

Our numerical result gives Itot — 1-085 (2) irL. Note that the value of Itot is linearly proportional to the width L and 
independent of boundary conditions. The numerical results of Gl (6) for both PJC and AJC with L = 20 are shown 
as circles in Fig. 3. There are two remarks on Fig. 3 worthy to be mentioned: (i) The differences in Gl (9) between 
AJC and PJC for L — 6 and 20 are shown in the up-left corner of Fig. 3, and it indicates the difference vanishes very 
quickly as L increases, (ii) By comparing with G^ (9) shown by the solid line, we know Gl (9) converges to Goo {&) 
very quickly as L increases. 

For the system with AJC, the critical point is one of the zeros for both even and odd L, and the cumulative 
distribution of zeros in the interval < 9 < #p=i im in is contributed solely by the zero mode p = 0. In this interval of 
9, we have Gl (9) = l p= o (9) /2l to t- By completing the integration of Eq. (13) for p = 0, we have 




(15) 

This result yield <?l (0) = 0.254(0) /L. Thus, by comparing Eq. (15) with Eq. (11) we may conclude the following 
properties for the system with AJC: (i) The system exhibits a first-order transition for finite L. (ii) The transition 
temperature always locates at the critical point (3 C of 2D Ising transition independent of the value of L. (Hi) The 
latent heat released by the transition is inversely proportional to the width L. 
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To obtain the exact amount of the latent heat released in the transition, we can measure the discontinuity of the 
internal energy at the transition point. The dimensionless internal energy density is defined as = dfi,/d(3. By 
rewriting the free energy density of Eq. (1) as 



h= ~\ In (4z) - -L £ [- to F (p, (3) + J(p, /?)] , (16) 

with 



2L 

p 



I(p, 0) = J* g In [1 - F(p, (3) cos 4>] , (17) 
F(p,0) = z(z 2 + l-cos^ , (18) 

we may obtain by first completing the integration of Eq. (17) and then performing the derivatives with respect to (3. 
However, to avoid any ill-defined result we have to be very cautious of continuity for the integrand. For conncreteness, 
we proceed first with the derivative and then the integration. By using the technique of contour integration, we have 

81 ( 1 \ 8 

Note that the above equality holds only when the condition 1 — F 2 (p,f3) > is met, and this condition ensures the 
continuity of the integrand. Then, we have 

p 

with 

P e (p,/?) = - — 1= = -jLbxF(p,0). (21) 
^1-F 2 (p, J) dp 

For PJC, the functional values of F(p, 0) satisfy the inequality < F(p, 0) < 1, and the internal energy of Eq. (20) 
is a continuous function of (3. For AJC, the functional values of F(p,(3) also satisfy the condition 1 — F 2 (p,f3) > 
with one exceptional point locating at p = and (3 = (3 C for which we have 1 — F 2 (0,f3 c ) = O^tJ- Thus, for AJC 
we may expect to have the mismatch of the left and right derivatives of the free energy at [3 C . The discontinuity 
A e defined as the differences between €l (f3 c — e) and ((3 C + e) exactly is the latent heat per site released by the 
transition. Since the zero mode of Eq. (21) is solely responsible for the discontinuity, we have 

A e = -L lim [D e (0, 13 C - e) - D e (0, 0c + e)] ■ (22) 

This yields A e = 4/L. Thus, the latent heat scales as 1/L, and this agrees with the result obtained from the analysis 
of the zeros. 

In summary, we analyze the distribution of Fisher zeros in the complex sinh (2rj) plane for an infinitely long Ising 
strip with finite width L joined periodically or antiperiodically. For periodically joined circumferece, the the leading 
finite-size scaling behavior for the distribution of Fisher zeros give the exact values of the correlation-length and shift 
exponents, v and A zero , as v = 1 and \ zero = 2. For antiperiodically joined circumference, the system contains 
the interface, and the fluctuation of the interface renders the critical point of Ising transition to be non-analytic 
in thermodynamic quantities for any width L. The nature of the non-analytic associated with the critical point is 
identified as the first-order phase transition for finite L. Because of the existence of the interface, the system can not 
be identified as an effective one-dimensional system with short-range Hamiltonian, and the appearence of the first- 
order transition does not contradict with the Mermin-Wigner theorem^|.The latent heat released by the transition 
scales exactly as 1/L with the the exact amount 4/L. Thus, as L increases, the strength of the first-order transition 
decreases, and eventually the system exhibits the second-order phase transition in the limit of L — > oo. 
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FIG. 1: The curves of 8 P versus cfi/n for various p- modes with L = 6 and periodic (solid lines) and antiperiodic (dot lines) 
boundary conditions. The number on a curve is the p-value 
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FIG. 2: The distributions of Fisher zeros in the complex z plane for infinitely long Ising strip with width L — 5 for (a) periodic 
and (b) antiperiodic boundary conditions and width L = 6 for (c) periodic and (d) antiperiodic boundary conditions. The 
opening angle is min = cos _1 [cos 2 (7r/2L)] 



FIG. 3: The cumulative distribution of Fisher zeros G (8) as a function of 9. The circles are the results for infinitely long Ising 
strip with width L = 20 for periodic and antiperiodic boundary conditions, and the solid line corresponds to the results from 
the integration of the zeros density given by Ref. [8]. The differences in Gl {&) between the antiperiodic and the periodic 
boundary conditions for L — 6 (dot line) and 20 (solid line) are shown in the up-left corner. 



